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Waveform  modeling  of  the  leleseismic  long-period  body-wave  phases  of  shallow  earthquakes  has  proven  quite 
effective  in  determining  source  parameters  for  events  of  magnitude  larger  than  six.  Unfortunately,  these  leleseismic 
phases  become  too  weak  for  smaller  events,  and  regional  data  must  be  used  which  are  generally  much  more 
complicated.  This  is  because  the  crust-mantle  system  acts  as  a  leaky  waveguide  and  a  large  number  of  rays  are  required 
to  model  the  observations.  In  this  report  we  review  a  procedure  for  the  systematic  determination  of  source  parameters 
from  regional  body  waves.  A  least-squares  inversion  technique  is  used  which  is  based  on  a  cross-correlation  of  the  data 
and  a  synthetic  seismogram.  The  synthetic  waveforms  are  a  linear  combination  of  those  for  the  three  fundamental  types 
of  fault.  A  set  of  profiles  of  the  synthetic  waveforms  is  presented.  Although  the  synthetic  waveforms  are  generated  for  a 
model  developed  for  the  western  United  States,  this  model  seems  to  be  adequate  for  most  continental  regions.  The 
inversion  is  parameterized  in  terms  of  fault  strike,  dip  and  rake;  these  parameters  are  relatively  insensitive  to  small 
changes  in  crustal  structure.  Several  examples  are  presented. 


1.  Introduction 

A  considerable  amount  of  effort  has  been  ex¬ 
pended  to  determine  the  source  parameters  of 
moderate-size  earthquakes,  although  such  de¬ 
terminations  can  be  beset  with  difficulties.  Ideally, 
a  large  amount  of  information  can  be  derived  from 
the  modeling  of  long-period  body  waves  (see 
Helmberger,  1974;  Langston  and  Helmberger, 
1975).  Unfortunately,  if  an  earthquake  is  too  small 
to  be  well  recorded  teleseismically,  which  is  the 
case  for  many  events  in  the  magnitude  range  be¬ 
tween  5  and  6,  the  fault-plane  orientation  must  be 
constrained  by  local  short-period  data  and  the 

*  Contribution  No.  3706,  Division  of  Geological  and  Planetary 
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seismic  moment  usually  cannot  be  determined  un¬ 
ambiguously.  The  World-Wide  Standard  Seismo¬ 
graph  Network  (WWSSN)  supplemented  by  other 
long-period  stations  and  arrays  provides  suffi¬ 
ciently  dense  coverage,  in  that  most  moderate-size 
earthquakes  occurring  in  continental  regions  will 
produce  some  on-scale  records  of  long-period  body 
waves  at  regional  distances.  In  the  regional  dis¬ 
tance  range  (1-12°),  the  waveguide  properties  of 
the  crust  produce  complicated  body-wave  signals; 
however,  in  most  cases  the  long-period  waveform 
is  quite  distinctive  and  sufficiently  insensitive  to 
details  of  crustal  structure  to  allow  the  separation 
of  the  source  and  structural  information. 

In  this  paper  we  review  a  procedure  for  extract¬ 
ing  the  source  parameters  of  moderate-size  earth¬ 
quakes  from  data  for  long-period  regional  phases. 
The  technique  involves  an  iterative  inversion  pro- 


0031-9201/82/0000-0000/S02.75  ©  1982  Elsevier  Scientific  Publishing  Company 

63  03  14  139 


186 


cess  which  minimizes  the  difference  between  a 
synthetic  seismogram  and  the  observation.  The 
synthetic  waveforms  are  constructed  using  Green's 
functions  computed  for  a  single,  very  simple  struc¬ 
ture.  These  Green’s  functions  appear  to  be  an 
adequate  model  for  most  continental  regions,  thus 
allowing  a  quick  and  approximate  determination 
of  fault  parameters.  The  inversion  is  parameterized 
in  terms  of  strike,  dip  and  rake.  The  number  of 
inversion  parameters  has  been  minimized  so  that 
inadequacies  in  the  Green’s  functions  are  not  over¬ 
emphasized.  Obviously,  the  structural  model  is 
more  appropriate  for  certain  regions  than  for 
others,  so  the  inversion  parameters  chosen  are 
those  which  are  most  robust.  The  main  advantage 
of  the  technique  is  that  it  requires  only  a  small 
data  set.  The  general  usefulness  of  this  technique  is 
illustrated  by  inverting  regional  data  for  earth¬ 
quakes  occurring  in  the  United  States,  northern 
Canada  and  southern  Europe. 


2.  The  Green's  functions 

The  techniques  for  constructing  the  Green’s 
functions  have  been  discussed  in  detail  elsewhere 
(Helmberger  and  Engen,  1980;  Wallace  et  al„ 
1981).  What  is  discussed  here  is  some  simplifying 
approximations  and  their  applicability.  The  wave¬ 
form  of  interest  is  that  part  of  a  regional  seismo¬ 
gram  which  arrives  before  the  S-wave.  This  wave¬ 
form  is  referred  to  as  P„,.  It  is  simplest  to  discuss 
P„,  in  terms  of  rays  which  are  in  a  waveguide.  The 
first  part  of  P„,  is  dominated  by  P  waves  (Pn)  and, 
moving  back  into  the  record,  the  waveform  con¬ 
tains  progressively  more  SV  (PL)  contributions. 
The  SV  energy  corresponds  to  rays  which  are 
reflected  within  the  crust  and  undergo  subsequent 
mode  changes  at  the  free  surface  and  the  Moho. 
The  interference  of  all  the  rays  gives  rise  to  the 
waveform;  parameters  such  as  the  crustal  thick¬ 
ness  or  the  velocity  contrast  between  the  crust  and 
the  Moho  control  the  waveform  dispersion.  Figure 
1  shows  some  typical  P„,  waveforms. 

The  Green’s  functions  are  constructed  by  sum¬ 
ming  generalized  rays  for  a  point  shear  disloca¬ 
tion.  As  an  example,  consider  the  following  equa¬ 
tion  for  the  vertical  displacement,  in  cylindrical 


Fig.  1.  The  horizontal  component  of  motion  for  four 
moderate-size  earthquakes.  The  magnitude  and  distance  to  the 
recording  station  are  given  for  each.  Clockwise  from  the  upper 
left:  (I)  September  12,  1966,  Truckee,  CA,  strike-slip  event;  (2) 
June  13,  1965,  southwest  Turkey,  normal  event;  (3)  August  1. 
1975,  Oroville,  CA.  normal  event;  and  (4)  December  10.  1967. 
off  the  coast  near  Cape  Mendocino,  CA.  strike-slip  event. 


coordinates 

w(r,z,0,  f )  =  ( A/0/4rrp)|  Z>(0  *  £ 

(1) 

where  D(t)  is  the  far-field  time  history,  p  the 
source  region  density,  and  M0  the  seismic  moment. 
The  summation  process  adds  the  contributions  of 
the  three  fundamental  types  of  fault;  W,  are  the 
Green’s  functions  for  vertical  strike-slip,  vertical 
dip-slip  and  45°  dip-slip  step  dislocations,  re¬ 
spectively.  At  are  coefficients  determined  by  the 
source  orientation  and  are  given  by 

,4,(0,  A,  5)  =  sin  20  cos  A  sin  5 

+  ( I  /2)cos  20  sin  A  sin  25 
/42(0,A,5)  =  cos  0  cos  A  cos  5  —  sin  0  sin  A  cos  25 


,4j(0,A.5)  =  (l /2)sin  A  sin  20  (2) 
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where  6  is  the  receiver  azimuth  from  the  end  of  the 
fault  plane,  \  is  the  rake  angle,  and  S  is  the  dip 
angle.  The  total  displacement  is  then  the  sum  of 
the  displacements  from  each  ray.  The  number  of 
rays  which  are  used  is  determined  by  the  structure 
and  the  pass  band  of  the  observation.  In  our 
previous  work  we  have  shown  that  a  single  layer, 
corresponding  to  the  crust,  over  a  half-space  man¬ 
tle  is  a  sufficient  structural  model  for  regional 
long-period  records  to  allow  the  extraction  of  the 
source  parameters  of  moderate-size  earthquakes. 

As  an  example.  Fig.  2  shows  a  comparison  of 
the  synthetic  waveforms  and  the  records  of  the 
1966  Truckee,  CA  earthquake  (which  will  be  dis¬ 
cussed  in  more  detail).  The  synthetic  waveforms 


were  constructed  from  Green’s  functions  com¬ 
puted  for  the  crustal  model  in  Table  I.  and  the 
fault  orientation  was  determined  by  inversion  of 
the  regional  data.  In  this  case  the  source  time 
function  and  moment  were  determined  by  other 
methods,  so  that  the  synthetic  amplitudes  can  be 
viewed  as  predictions.  The  numbers  on  the  traces 
are  maximum  peak-to-peak  amplitudes.  The  only 
noticeable  difference  between  the  data  and  the 
synthetic  waveforms  concerns  the  high-frequency 
content,  which  could  be  caused  by  several  factors: 

( 1 )  the  effects  of  attenuation  within  the  crust  have 
not  been  added  to  the  synthetic  waveforms,  and 

(2)  the  very  sharp  boundaries  in  our  model  are 
efficient  in  trapping  short-period  energy.  Overall, 


0 


Fig.  2.  Vertical  p„,  waveforms  from  the  Truckee  earthquake.  The  star  denotes  the  epicenter.  The  data  are  the  top  trace  at  each  station, 
the  trace  below  is  the  synthetic  fit.  The  strike-slip  mechanism  has  two  nodal  planes  which  project  through  the  stations  TUC  and  BOZ. 
To  the  right  of  each  trace  is  the  observed  or  predicted  amplitude  (on  the  basis  of  a  moment  of  0.8x  102'  dyne  cm)  in  10  ’  cm. 
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TABLE  I 
Crustal  model 


<v 

v. 

Density 

Layer 

(km  s  ') 

(km  s') 

(g  cm-3) 

thickness 

(km) 

6.2 

3.5 

2.7 

32 

8.2 

4.5 

3.4 

the  fit  of  the  synthetic  waveforms  to  the  data 
justifies  the  use  of  the  simple  model,  and  the 
high-frequency  content  of  the  synthetic  waveforms 
does  not  affect  our  ability  to  determine  the  source 
parameters. 

The  only  expensive  or  complicated  part  of  this 
modeling  process  is  the  actual  generation  of  the 
Green’s  functions.  For  this  reason.  Figs.  3  and  4 
are  presented,  which  give  profiles  of  the  vertical 
and  radial  responses  for  the  three  fundamental 
faults.  Most  earthquakes  which  produce  an  on- 
scale  P„,  at  long-period  WWSSN  stations  have 
similar  time  functions  (which  is  simply  a  reflection 
of  event  size).  In  Figs.  3  and  4  a  trapezoid  with  a 
1  s  rise,  a  1  s  top  and  a  1  s  fall  is  used  for  the 
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Fig.  3.  Theoretical  displacement  profiles  for  the  vertical  compo¬ 
nent.  The  Green's  functions  were  computed  from  the  model 
presented  in  Table  I  and  have  been  convolved  with  a  source 
time  function  represented  by  a  trapezoid  (r,  =  1. 12  =  1. 1,  =  I), 
a  triangular  filter  (2  s  rise  and  fall),  and  the  response  function 
for  a  WWSSN  long-period  instrument. 


RSS 

RD5 

R' 

•  '  4.5  v 

4.5  • 

5.4  . 

54 

deq 

*•  •  1 

-  y;  63.- 

6.3  •  ‘  ' 

7  2  . 

72 

-  '  ■•••-  8.'  - 

8i  -  -  • 

ao-. 

90  - ;  * 

•  .  9.9  . 

A  *\  " 

9.9  -  • 

10.8  •■■■-  ' 

’  -cV'ii.r-'!' 

; . 

1(7  - 

*  1 

1 

Fig.  4.  Theoretical  displacement  profiles  for  the  radial  compo¬ 
nent.  The  Green's  functions  are  computed  every  100  km.  They 
have  been  convolved  with  the  time  function,  instrument  re¬ 
sponse  and  filter  described  in  Fig.  3. 


far-field  time  history.  The  displacements  have  also 
been  convolved  with  the  response  function  for  a 
15-100  instrument.  Because  of  the  differences  in 
high-frequency  content  between  the  data  and  the 
synthetic  waveforms,  the  displacement  responses 
have  also  been  filtered.  The  filter  has  an  impulse 
response  represented  by  a  triangle  which  has  a  2  s 
rise  and  fall.  When  comparing  these  displacements 
with  the  data,  the  observations  should  be  similarly 
filtered.  Once  the  response  for  the  three  funda¬ 
mental  faults  is  known  any  seismogram  can  be 
constructed  by  a  linear  combination  of  them. 

The  displacements  in  Figs.  3  and  4  were  com¬ 
puted  for  a  source  depth  of  8  km.  Varying  the 
source  depth  between  5  and  15  km  has  only  a 
small  effect  on  the  waveform.  This  is  easily  under¬ 
stood  by  considering  that  to  first  order,  a  change 
in  source  depth  affects  only  the  travel  time  of  the 
first  segment  of  any  ray.  Figure  5  shows  a  com¬ 
parison  of  the  synthetic  waveforms  at  1000  km  for 
three  different  source  depths.  After  doubling  the 
source  depth  (from  8  to  16  km)  the  essential 
character  of  the  waveform  is  still  preserved  and 
the  source  information  is  retrievable.  In  contrast,  a 
similar  change  in  crustal  thickness  would  affect  the 
travel  time  of  each  leg  of  a  given  ray,  hence 
significantly  changing  the  waveform  dispersion 
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Fig.  5.  P„,  waveforms  al  1000  km  for  three  different  source 
depths. 


(Wallace  and  Helmberger,  1980).  The  insensitivity 
of  the  displacements  to  source  depth  allows  the 
responses  in  Figs.  3  and  4  to  be  used,  at  least  in  a 
qualitative  fashion,  to  determine  the  source 
parameters  of  most  crustal  earthquakes. 

The  only  other  major  question  concerning  the 
applicability  of  the  displacements  presented  relaies 
to  the  structure  used  in  their  calculation.  Experi¬ 
ence  indicates  that  the  simplistic  model  adopted  is 
justified.  The  model  in  Table  I  is  an  average  devel¬ 
oped  for  the  western  United  States,  although  it 
appears  to  be  sufficient  for  most  continental  re¬ 
gions  in  the  world.  The  waveform  dispersion  is 
dependent  on  crustal  thickness  and  on  the  contrast 
between  the  upper-mantle  P  velocity  and  the  mean 
crustal  P  velocity,  so  obviously  for  regions  having 
anomalous  crustal  structure,  such  as  the  Tibetan 
Plateau,  the  responses  in  Figs.  3  and  4  would  be 
inadequate.  Also,  the  use  of  the  half-space  to 
approximate  the  upper  mantle  must  cease  to  be 
valid  at  some  point;  at  some  distance  a  significant 
amount  of  energy  will  be  present  in  the  form  of 
diving  rays  which  have  turned  in  the  mantle.  For 
most  continental  regions  this  distance  appears  to 
be  —  12°.  If  diving  rays  are  present,  the  ratio  of 
the  Pn  to  PL  amplitudes  should  differ  from  the 
synthetic  value  (Wallace  and  Helmberger,  1981). 
Also,  the  predicted  amplitude  should  begin  to 
diverge  from  that  observed,  owing  to  attenuation 
in  the  upper  mantle.  Figure  6  shows  a  profile  of 
the  Truckee  earthquake  records;  the  amplitudes 
have  been  corrected  for  the  azimuthal  radiation 
pattern  and  the  polarities  have  been  adjusted  to 
show  a  smooth  dispersion  pattern.  To  the  right  is  a 
profile  of  displacements  for  a  strike-slip  earth¬ 
quake  (such  a  profile  can  be  constructed  from 
Fig.  3).  Note  that  there  does  not  appear  to  be  a 


Fig.  6.  Truckee  earthquake  waveforms  corrected  for  horizontal 
radiation  pattern  and  plotted  as  a  function  of  distance.  The 
maximum  amplitude  is  shown  to  the  right  of  each  trace.  Note 
that  the  stations  BOZ  and  TUC  are  very  close  to  nodes. 

systematic  break-down  in  waveform  shape  or  am¬ 
plitude  over  the  distance  range  4-12°.  Also  note 
that  the  stations  BOZ  and  TUC  are  nearly  nodal, 
and  thus  their  recordings  of  amplitude  are  not 
particularly  reliable. 


3.  Inversion  technique 

The  ability  to  determine  the  source  parameters 
of  an  earthquake  by  comparing  an  observed  with  a 
predicted  waveform  depends  on  the  assessment  of 
the  quality  of  fit.  In  a  previous  paper  (Wallace  et 
al.,  1981)  we  presented  a  least-squares  waveform- 
inversion  technique  which  makes  use  of  an  error 
function  determined  by  the  cross-correlation  of  a 
long-period  seismogram  and  a  synthetic  waveform; 

Vf  (3) 

where  /  is  the  observed  record,  g  is  the  synthetic 
waveform,  and  the  integral  is  a  zero-lag  cross-cor¬ 
relation.  The  limits  of  integration  are  the  time 
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length  of  the  window  in  which  the  waveforms  are 
correlated.  The  denominator  serves  to  normalize 
both  the  data  and  the  synthetic  waveform.  This 
normalization  process  renders  the  error  function 
insensitive  to  the  absolute  amplitudes.  To  mini¬ 
mize  the  error,  which  corresponds  to  maximizing 
the  correlation.  /  and  g  are  allowed  to  optimally 
align  themselves  with  regard  to  waveform;  /  and  g 
are  aligned  a  priori  in  time  by  matching  first 
breaks  and  ignoring  absolute  travel  time.  The  error 
function  can  be  rewritten  by  considering  that  the 
synthetic  seismogram  can  be  constructed  using 
waveforms  for  the  three  fundamental  faults.  In 
this  case  there  is  a  summation  of  the  cross-correla¬ 
tions  between  the  observed  and  each  of  the  funda¬ 
mental  faults.  For  a  given  range,  the  cross-correla¬ 
tions  are  constant  and  errors  can  be  minimized  by 
varying  the  4,'s  in  eq.  (2).  In  other  words,  once  the 
cross-correlations  are  computed  the  source  orien¬ 
tation  is  determined  iteratively  and  only  the  con¬ 
stants  have  to  be  recalculated.  For  a  detailed  dis¬ 
cussion  of  the  inversion  technique  see  Wallace  et 
al.  (1981). 

Once  the  source  orientation  is  fixed,  the  mo¬ 
ment  of  an  earthquake  can  be  determined  by 
comparing  the  amplitudes  of  the  synthetic 
waveform  and  the  observations.  Adopting  the  units 
of  Helmberger  and  Malone  (1975),  and  expressing 
the  range  in  km,  time  in  s,  density  in  g  cm  \ 
velocity  in  km  s'1,  the  moment  in  dyne  cm,  and 
displacement  in  cm.  yields 


M,  =  4.P  X  1<W  -df  ar^Pllt^-  )  (4) 

synthetic  amplitude  / 


A  moment  can  be  determined  by  comparing  the 
maximum  peak-to-peak  amplitudes  for  any  time 
window  used  for  the  correlation.  It  has  been  found 
that  a  moment  should  be  determined  for  a  few 
peaks  at  a  given  station.  The  ratio  of  the  moment 
at  each  station  to  the  mean  is  a  measure  of  the 
amplitude  stability.  In  general,  the  moments  de¬ 
termined  from  are  in  very  good  agreement  with 
those  determined  teleseismically,  using  an  assump¬ 
tion  of  f*  =  1  s. 


4.  Examples 

We  have  inverted  the  P„,  waveforms  from  five 
earthquakes  to  demonstrate  the  utility  of  the 
technique.  Three  of  the  earthquakes  occurred  in 
the  western  United  States.  Two  other  earthquakes, 
one  in  Baffin  Bay  in  the  Artie  and  the  other  in 
Turkey,  have  been  included  to  demonstrate  that 
the  Green's  functions  are  not  unique  to  the  United 
States.  Both  dip-slip  and  strike-slip  mechanisms 
are  represented  in  the  suite  of  examples. 

4.1.  Truckee,  CA  (12/9/67) 

The  Truckee  earthquake  was  a  strike-slip  event 
at  10  km  depth  which  produced  excellent  regional 
records  but  very  few  teleseismic  body-wave  re¬ 
cords,  as  is  typical  of  moderate-size  strike-slip 
events.  The  Truckee  earthquake  tmh  =  5.7)  has 
been  studied  by  numerous  authors  (Ryall  et  al.. 
1968;  Tsai  and  Aki,  1970;  Burdick.  1977).  making 
it  a  good  test  case.  Tsai  and  Aki  (1970).  from 
first-motion  studies  and  modeling  of  the  surface 
waves,  determined  this  event  to  be  pure  strike- slip 
on  a  fault  plane  striking  N44°E  and  dipping  80°SE. 
The  surface  wave  moment  was  determined  to  be 
0.83  X  1025  dyne  cm.  Figure  2  shows  the  locations 
of  the  epicenter  and  the  recording  stations,  and  the 
filtered  data  for  Truckee.  Also  shown  are  the 
synthetic  waveforms  determined  from  the  inver¬ 
sion  results.  Note  that  the  stations  BOZ  and  TUC 
are  very  nearly  nodal.  The  inversion  yields  a  mech¬ 
anism  which  is  very  similar  to  Tsai  and  Aki's 
(1970);  a  strike  of  N43°E,  a  dip  of  76°SE,  and  a 
rake  of  —11°.  The  only  significant  difference  is 
the  slight  dip-slip  component  in  our  solution, 
which  is  also  acceptable  on  the  basis  of  the  first- 
motion  data.  The  moment  determined  from  the  P„, 
waveforms  is  0.87  X  1025  dyne  cm.  which  is  in 
excellent  agreement  with  the  result  obtained  by 
Tsai  and  Aki  (1970). 

4.2.  El  Golfo,  Mexico  (7/8/66) 

The  El  Golfo  earthquake  (mb-  6.3,  A/s  =  6.3) 
was  a  strike-slip  event  which  occurred  near  the 
mouth  of  the  Colorado  River  at  the  northern  end 
of  the  Gulf  of  California.  Ebel  et  al.  (1978)  de- 
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4.3.  Oroville.  CA  (1/8/75) 


termined  the  fault  plane  to  be  striking  E140°S, 
dipping  85°  to  the  southwest,  and  with  a  rake  of 
183°,  and  determined  the  depth  of  the  event  to  be 
10  km.  Using  teleseismic  long-period  P  waves  they 
determined  a  moment  of  5.0  X  1025  dyne  cm. 

El  Golfo  is  about  the  maximum-size  event  which 
can  be  used  in  the  inversion  technique.  The  P„/ 
records  are  barely  on-scale  at  the  stations  used. 
Figure  7  shows  the  locations  of  the  epicenter  and 
the  recording  stations,  and  the  waveforms.  In  this 
case  the  time  function  is  a  triangle  with  a  2  s  rise 
and  fall;  the  long-period  impulse  is  a  reflection  of 
the  event  size.  The  long-period  time  function  al¬ 
lows  the  use  of  a  filter  to  be  dispensed  with.  In 
Fig.  7,  shown  below  the  observed  waveforms  are 
the  synthetic  waveforms  for  the  inversion  solution, 
namely,  for  a  strike  of  E137°S,  a  dip  of  81°,  and  a 
rake  of  175°.  The  inversion  solution  is  in  good 
agreement  with  the  results  of  Ebel  et  al.  (1978). 
The  moment  determined  from  the  Pn/  waveforms  is 
4.6  X  1025  dyne  cm. 


Fig.  7.  Data  and  synthetic  waveforms  for  the  El  Golfo  earth¬ 
quake.  The  map  gives  the  location  of  the  event  (star)  and  the 
recording  stations.  Along  each  trace  is  the  ratio  of  the  station 
moment  to  the  average  moment. 


The  Oroville  earthquake  (M  =  5.6)  was  a  nor¬ 
mal  faulting  event  and  is  interesting  because  the 
surface-wave  (Hart  et  al„  1977)  and  body-wave 
(Langston  and  Butler.  1976)  analyses  yield  sub¬ 
stantially  different  moments.  Langston  and  Butler 
(1976)  determined  the  strike  to  be  180°,  with  a  dip 
of  65°  and  a  rake  of  —70°.  Their  moment  de¬ 
termination  is  5.7  X  1024  dyne  cm.  Hart  et  al. 
(1977)  suggested  that  the  surface  waves  are  con¬ 
sistent  with  this  body-wave  mechanism  but  that 
the  moment  is  larger  by  a  factor  of  3  (1.9x  102* 
dyne  cm).  Figure  8  shows  the  locations  of  the 
event,  of  the  stations  used  in  the  inversion  analy¬ 
sis,  and  the  filtered  data  and  synthetic  waveforms. 
The  inversion  solution  has  shifted  the  mechanism 
to  a  strike  of  204°.  dipping  66°  with  a  rake  of 
—  85°.  This  new  solution  violates  only  a  few  first 
motions,  but  the  aftershock  trend  tends  to  support 
the  180°  strike.  The  moment  determined  from  P,(/ 
data  is  6.9  X  1024  dyne  cm. 

4.4.  Baffin  Bay.  Canada  (4/9/63) 

The  Baffin  Bay  earthquake  ( M  —  5.9)  was  a 
normal  event  associated  with  a  continental  margin 
area.  The  travel  path  to  each  of  the  stations  used 
in  the  inversion  includes  portions  of  crustal  and 
oceanic  regions,  which  makes  it  an  ideal  event  to 
test  the  applicability  of  the  Green’s  functions.  Liu 
and  Kanamori  (1980)  modeled  the  body  waves 
and  determined  a  fault-plane  solution  with  a  strike 
of  98°,  a  dip  of  66°N,  and  a  rake  of  —  103°.  The 
location  of  the  event,  and  the  filtered  data  and 
synthetic  waveforms  are  shown  in  Fig.  9.  The  in¬ 
version  solution  has  a  mechanism  striking  74°. 
dipping  66°.  and  with  a  rake  of  —  100°.  Again  the 
only  parameter  which  is  appreciably  different  from 
the  results  of  teleseismic  analysis  is  the  strike.  In 
any  case  the  inversion  solution  is  acceptable,  and 
the  difference  of  20°  in  strike  can  be  considered  as 
the  resolution  for  dip-slip  events. 

4.5.  Turkey  (13/6/65) 

The  Turkish  event  ( A/  =5.1)  was  a  shallow, 
normal  event  which  occurred  in  southwest  Turkey 
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Fig.  8.  Filtered  data  and  synthetic  waveforms  for  the  Oroville  earthquake.  At  all  the  stations  except  GOL  both  the  vertical  (the  first 
trace  pair)  and  radial  components  are  shown. 


in  a  region  of  north-south  extension.  McKenzie 
(1972)  used  first-motion  data  to  determine  a  pure 
normal  mechanism  with  a  strike  of  101°  and  dip¬ 
ping  70°  to  the  south,  although  it  is  not  well 
constrained.  There  were  three  WWSSN  stations  at 
regional  distances  which  could  be  used  in  the 
inversion  process.  Figure  10  shows  the  locations  of 
the  event  and  the  recording  stations.  The  filtered 


data  and  the  fit  of  the  synthetic  waveforms  are 
also  shown.  The  inversion  solution  (strike  131°. 
dip  68°,  rake  —88°)  is  consistent  with  the  first- 
motion  data,  although  it  differs  in  strike  from 
McKenzie’s  solution.  Again  the  three-station  solu¬ 
tion  is  quite  acceptable,  considering  the  quality  of 
the  first-motion  data. 
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BAfTtN  BAY  9/4/63 


Fig.  9.  Location  of  the  Baffin  Bay  earthquake  (star)  and  the  recording  stations.  Filtered  data  and  synthetic  waveforms  for  both  the 
vertical  and  radial  components  are  shown. 


Fig.  10.  Location  of  the  Turkey  earthquake  (star)  and  the  recording  stations.  Filtered  data  and  synthetic  waveforms  for  both  the 
vertical  and  radial  components  are  shown. 


5.  Discussion 

Determining  the  fault-plane  orientation  of 
moderate-size  earthquakes  is  often  a  frustrating 
experience,  owing  to  the  paucity  of  high-quality 
data.  Earthquakes  in  the  magnitude  range  S-6  are 
quite  important  and  are  often  the  only  “measura¬ 


ble”  expression  of  the  present  tectonic  environ¬ 
ment.  All  of  the  available  data  must  be  used  to 
extract  the  source  parameters  of  these  moderate- 
size  events,  and  the  modeling  of  Pnl  waveforms  can 
provide  a  valuable  constraint  in  this  process.  Every 
situation  will  probably  be  unique  and  it  is  difficult 
to  predict  which  data  set  will  be  the  most  defini- 
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tive.  Nevertheless,  it  appears  worthwhile  to  con¬ 
sider  the  P,w  inversion  separately  and  its  resolva¬ 
bility  dependence  on  source  orientation. 

The  inversion  process  discussed  requires  only  a 
small  data  set.  With  ideal  azimuthal  separation,  a 
data  set  comprisingjust  three  stations  (vertical  and 
radial  components)  can  yield  good  solutions  (it  is 
possible  to  construct  a  case  where  the  inversion  is 
unstable,  but  in  practice  this  has  never  happened). 
In  almost  all  cases  four  recording  stations  are 
sufficient.  P„,  should  be  polarized  in  the  vertical 
and  radial  planes.  Rotation  of  the  horizontal  com¬ 
ponents  for  a  number  of  events  indicates  that 
there  is  very  little  energy  on  the  tangential  compo¬ 
nent.  implying  little  contamination  from  such  ef¬ 
fects  as  multipathing.  Therefore,  it  is  usually  suffi¬ 
cient  to  take  the  largest  of  the  horizontal  compo¬ 
nents  to  be  the  radial  waveform  in  the  inversion. 
The  resolving  ability  of  the  inversion  (or  con¬ 
versely.  the  error)  depends  on  the  type  of  earth¬ 
quake.  The  experience  gained  by  considering  the 
examples  presented  in  the  last  Section  indicates 
that  the  mechanisms  of  strike-slip  earthquakes 
can  be  determined  quite  well  using  data  from 
relatively  few  stations.  The  strike  is  usually  de¬ 
termined  to  within  5°  of  that  determined  by  other 
methods.  The  rake  is  the  least  resolvable  parame¬ 
ter  for  strike-slip  events  and  can  vary  by  up  to 
15°  from  that  determined  by  first-motion  studies. 
The  mechanisms  of  dip-slip  earthquakes  are  more 
difficult  to  determine.  Although  the  dip  and  rake 
are  usually  determined  in  good  agreement  with 
other  studies,  the  strike  may  vary  by  up  to  20°. 
This  feature  is  illustrated  by  considering  a  45° 
dipping  normal  fault.  In  this  case,  most  regional 
stations  lie  within  the  compressional  region  of  the 
focal  sphere  and  any  given  azimuth  will  produce 
remarkably  similar  waveforms.  Fortunately,  dip- 
slip  events  have  rather  strong  teleseismic  P  waves, 
but  again  the  waveforms  are  all  the  same,  with 
little  dependence  on  azimuth.  In  this  case  the 
stations  lie  in  the  center  of  the  focal  sphere  and 
are  all  dilatational.  As  an  example  consider  Fig. 
11,  which  displays  the  focal  sphere  and  the  tele- 
seismic  waveforms  for  the  Oroville  earthquake. 
Note  that  the  P  waveforms  are  all  similar,  and  the 
synthetic  waveforms  (Langston  and  Butler.  1976) 
do  not  add  much  insight  into  determining  the 


each  trace  is  the  moment  (  ♦  IU:'  dvne  cmi  determined  at  that 
station.  Note  the  scatter  over  a  factor  of  .'  (figure  from  I  angx- 
lon  and  Butler.  1976). 

strike  of  the  fault.  Figure  12  shows  the  filtered 
regional  data  and  synthetic  waveforms  computed 
from  the  teleseismic  fault  parameter  determina¬ 
tions.  A  comparison  of  Figs  12  and  S  shows  that 
the  regional  data-inversion  solution  improves  the 
fit  of  the  synthetic  waveforms,  in  particular  PAS 
This  suggests  that  a  logical  approach  would  be  to 
invert  some  teleseismic  data  and  the  regional  data 
simultaneously.  Since  the  inversion  technique  re¬ 
lies  on  the  cross-correlation  of  data  and  a  s\n- 
thetic  waveform,  the  joint  inversion  is  quite  tracta¬ 
ble. 

The  higher  resolution  of  the  strike  slip  events  is 
actually  fortuitous.  Moderate-size  strike  slip 
earthquakes  rarely  produce  usable  teleseismic  P 
waves,  owing  to  their  inefficiency  in  radiating  en¬ 
ergy  straight  down.  On  the  other  hand,  strike  slip 
events  produce  very  good  regional  waveforms  This 
allows  the  inclusion  of  a  larger  data  set  in  the 
inversion,  and  hence  the  resolution  problem  is  at 
least  partly  resolved. 

It  is  reasonable  to  consider  what  effect  the 
structural  model  has  on  the  inversion  results  As  a 
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X  = -70° 


Fig.  12.  Filtered  P„,  waveforms  from  the  Oroville  earthquake.  The  synthetic  waveforms  were  computed  with  the  teleseisntic  fault-plane 
solution. 


test  of  the  insensitivity  of  the  fault  orientation  to 
small  changes  in  crustal  parameters  the  El  Golfo 
earthquake  was  reinverted  with  a  different  struc¬ 
ture.  The  crustal  thickness  was  reduced  to  24  km, 
the  source  depth  was  moved  to  12  km,  and  the  Pn 
velocity  reduced  to  7.8  km  s  '.  Although  the 
quality  of  the  fit  decreases  significantly  the  mecha¬ 
nism  returned  by  the  inversion  is  similar;  strike 


138°.  dip  82°,  and  rake  181°,  The  moment  in¬ 
creases  to  6.9  X  10 25  dyne  cm. 


6.  Conclusion 

It  is  possible  to  extract  the  source  parameters  of 
moderate-size  earthquakes  from  long-period  re- 


1% 


gionai  body  waves.  The  procedure  requires  the 
comparison  of  the  observed  with  a  synthetic  wave¬ 
form:  the  synthetic  waveforms  can  be  generated 
by  a  linear  combination  of  the  waveforms  of  the 
three  fundamental  faults,  shown  in  Figs.  3  and  4. 
Although  these  synthetic  waveforms  are  for  a  sim¬ 
ple  model  the  inversion  parameters  (fault  strike, 
dip  and  rake)  are  fairly  insensitive  to  small  changes 
in  crustal  thickness.  Pn  velocity  and  mean  crustal 
velocity.  This  allows  this  single  set  of  Green’s 
functions  to  be  used  for  most  continental  earth¬ 
quakes.  The  inversion  procedure  requires  only  a 
small  data  set,  and  is  particularly  ideal  for  strike- 
slip  earthquakes. 
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